Testing the neutrality of matter by acoustic means in a spherical resonator 
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I. INTRODUCTION 

The most recent experimental value of the elementary 
charge e is p] 

e = 1.602 176 487(40) x 10~ 19 C 

where 0.000 000 040 x 10~ 19 C is the standard uncertainty, 
equivalent to 2.5 x 10 -8 relative standard uncertainty. It 
is commonly accepted that the charges of the electron 
and of the proton have equal magnitude, and this fact is 
usually associated with the principles of charge conser- 
vation and of baryon and lepton conservation, as well as 
with the symmetry between particles and antiparticles. 
However, as the charge symmetry lacks a corresponding 
conservation law [2], there is no compelling requirement 
for its validity: it rests uniquely on experimental bases. 
Nevertheless, and quite obviously, the structure of the 
elementary charges is of paramount importance for the- 
oretical physics since 1931, when Dirac pointed out that 
the existence of magnetic monopoles would imply charge 
quantization [3J . The uniqueness of the values of the ele- 
mentary charges for both leptons and baryons implies a 
novel connection, beyond Standard Model, between the 
two families of particles, suggesting, for example, the ex- 
istence of hypothetical lepto-quark particles which decay 
in both leptons and quarks [3]. 
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In the last century, several measurements have been 
performed in order to test the symmetry of the elemen- 
tary charges, using different laboratory approaches; to 
this purpose, most experiments probe the neutrality of 
common matter. The first laboratory method used a gas 
flowing out of a container whose electrostatic potential 
could be measured [5H7]. These measurements were ini- 
tially spurred by Einstein's suggestion that a small differ- 
ence between electron and proton charges could account 
for the magnetic field of the Earth [5] ; it was noted also 
that a small deviation from the neutrality of matter could 
explain the expansion of the universe [5J [5] ; these phe- 
nomena could also be due to a nonzero neutron charge 
6J. By now, both hypotheses have been ruled out by 
experiment. In the beam deflection method, the electric 
properties of a single species are tested: measurements 
are performed on molecules [H ITMT3"] and neutrons [T4T 
IT7] . A modern version of Millikan's experiment applies 
an electrostatic force to magnetically suspended small 
bodies [TS]; this method was initially devised to detect 
the presence of free quarks in common matter, but it 
allows also to put limits on the neutrality of matter |T9j- 
I2T] . The latest laboratory method applies an alternating 
electric field to a gas contained inside an acoustic cav- 
ity; if the gas molecules carry an electric charge, a sound 
wave is generated. This method is an original idea by 
Dylla and King [22] , and is often cited as one of the best 
results on this topic [23] • However, in the following we 
show that the results of the experiment were incorrect, as 
they were affected by a few errors. Among the planned 
experiments, we want to cite the proposed use of a tor- 
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sion balance [21] , and of atom interferometry [55] . Limits 
on charge asymmetry are cast also by model-dependent 
astrophysical methods [2DJ. Moreover, one should men- 
tion a test on the charges of positrons and antiprotons, 
derived from measurements of their cyclotron resonance 
frequencies and from spectroscopic data }27j . Limits have 
been set also on the neutrality of the neutrino by astro- 
physical methods |28j . 

The results of all the experiments performed so far 
are consistent with the usual view that matter is neu- 
tral, and have pushed further and further the limits on 
charge asymmetry between proton and electron and on 
the charge of the neutron. The results are best described 
in terms of two parameters 

% + <7c , <7n 

e p - e = and e n = — , 

e e 

where q v , q e and q n are the electric charges of protons, 
electrons and neutrons, respectively. If a body containing 
Z protons, Z electrons and N neutrons is measured to 
be neutral in an experiment capable of a sensitivity Sq, 
one has 

\Ze p _ c + Ne n \e < Sq. 

In order to disentangle the two contributions, one has 
to perform independent measurements on at least two 
systems with different Z and/or N. A common approach 
assumes instead charge conservation in the /3 decay of 
the neutron n — > p + e~ + v (and the neutrality of the 
antineutrino); with these hypotheses e p _ = e n = e q and 

I I < S 1 ^ S 1 m p 
{eql ~ (Z + N)e ~ e m ' 

where m p and m are respectively the masses of the proton 
and of the molecule or the sample. 

The actual experimental situation is summarized in 
Tab. |TJ The measurements performed with the gas efflux 
method, leaving aside Ref. [3D], give a limit on e q of the 
order of 10 -21 , a value obtained from careful subtraction 
of spurious effects of uncertain origin. The beam deflec- 
tion method has, for molecules, a sensitivity one order of 
magnitude worst than the previous one; it heavily relies 
on a modelization of the beam profile and energetics and 
suffers from electric field spatial inhomogeneities, which 
exert forces on neutral molecules; it is difficult to think 
to improve this experimental set-up. Also the suspended 
body method has reached a sensitivity of the order of 
10~ 21 , which might be improved by a factor 10 by in- 
creasing the statistics [3T]; gaining a larger factor relies 
on the possibility of using bigger bodies, which is not triv- 
ial |21j . In order to try to improve the limits on charge 
asymmetry, we chose the experimental set-up of a gas 
filled electric capacitor used as an acoustic resonant cav- 
ity [22], which seemed less sensitive to systematic arte- 
facts and at the same time most passible of upgrades. We 
considered, in particular, that a better sensitivity could 
be reached i) by increasing the gas pressure and hence the 
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Table I. One sigma upper limits on the neutrality of matter 
from direct laboratory measurements. Methods: 1, gas efflux; 
2, beam deflection; 3, suspended body; 4, acoustic. All the 
experiments are consistent with the assumption that matter is 
neutral. The present work questions the quoted limit obtained 
by method 4. 



quality factor Q of the cavity, and ii) by improving the 
detection method. In fact, in the original experiment the 
pressure in the cavity was limited to 1 atm; moreover we 
judged that we would have benefitted by the use of more 
modern electronics. But there were also other reasons for 
repeating the 1973 experiment by Dylla and King: in the 
first place the authors did not realize that their resonant 
mode was almost frequency degenerate with another one, 
a situation that calls for a complicated analysis. More- 
over, we show that their calculation of the neutrality sig- 
nal was wrong by a factor about sixty, while the value of 
the calibration signal they quote was wrong by a factor 
about three. As a consequence, their limit on electron- 
proton charge difference is, in reality, not better than 
about 1CT 19 . 

In Sect. II we develop a simple mathematical descrip- 
tion of the electrically excited acoustic effects in a spheri- 
cal resonator; Sect. Ill is a detailed account of the exper- 
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imental set-up and of the method; Sect. IV-A describes 
the results of the polarization measurements used to cal- 
ibrate the apparatus; neutrality results are discussed in 
Sect. IV-B; conclusions follow. 



II. PRINCIPLE OF MEASUREMENT 




Figure 1. Scheme of a spherical acoustic resonator. Rigid 
walls at radiuses n and T2 enclose an oscillating fluid. 



To measure the neutrality of matter with the acous- 
tic method [22 , an oscillating voltage is applied to a gas 
filled spherical capacitor with internal and external ra- 
diuses r\ and r2, respectively - see Fig. [I] Forces are 
exerted on the molecules in the first instance through 
electric polarization, but also if matter is not neutral. In 
both cases, acoustic waves are generated, whose ampli- 
tude is measured by means of a sensor placed in gas at 
the external surface. In the following we show that the 
two mechanisms can be discriminated; hence the method 
is able to cast limits on the non neutrality of matter. 

A detailed analysis of the paper by Dylla and King [22] 
revealed inconsistencies in their Eqs. (2) and (3). We 
have independently solved the differential equation for 
the acoustic wave both analytically |36j and numerically. 
The most general analytical solution which takes into full 
account the dissipative forces is too complicated to be 
included in the present work. Only a simplified version 
is presented below, and only for the neutrality problem. 
Anyway, all the results obtained from the complete, the 
simplified and the numerical solutions are found to fully 
coincide. 



A. Electric forces in gases 

The force per unit volume on an homogeneous dielec- 
tric with electric permittivity e = e e r and a charge den- 
sity p subject to an electric field E is given by [3T] 



with p m the mass density. The second term is the 
force associated with electric polarization. The dielectric 



constant relates to mass density through the Clausius- 
Mossotti law 



1 



.4 



1 



3A 
I- A 



with 



A = 



3e ' 



where a is the electric polarizability and n = NaPtu/M 
is the number of molecules per unit volume, with M the 
molecular weight. This leads to 



dt 

Op m 



(e r -l)(e r + 2). 



If the gas has low density and bears a charge asymmetry, 
then 

Pm 



f = {Ze, 



iVe n 



Zm u + Nrrin 



eE- 



N A p r , 



aVE 2 . (1) 

2 M K ' 

The two terms can be distinguished from their depen- 
dence on the electric field: for example, if the electric 
field is modulated at a frequency v, the first term is as- 
sociated with the response of the gas at frequency v and 
allows the study of the charge asymmetry, while the sec- 
ond term generates a signal at frequency 1v which can 
be used for calibration of the apparatus. 



B. Sound wave equation 

To first order in the vibration amplitude, and neglect- 
ing viscosity, the equation of motion of a small portion 
of fluid is 



VP + f = Pm0 



dv 



(2) 



where P is the pressure, p m o the unperturbed mass den- 
sity and v the fluid velocity. In the same approximation, 
conservation of mass is written as 



dp„i 
dt 



(3) 



Since acoustic waves propagate with little or no heat ex- 
change between adjacent regions of the fluid, the left- 
hand side of this equation can be transformed to 



dp m _ ( dp m 
dt ~ \ dP 



dP , dP 



where the subscript S indicates an isoentropic process 
and ks is the adiabatic compressibility modulus of the 
fluid. By taking the divergence of Eq. ^ and the time 
derivative of Eq. ([3]) and summing, we obtain an inho- 
mogeneous scalar wave equation for the pressure P: 



V 2 P 



1 d 2 P 



V/ 



(4) 
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with wave speed c given by 
dP 



1 



ks PmO 



If the fluid is an ideal gas, then the Laplace's equation 
holds 



-yPo 

PmO 



jRT 
M ' 



(5) 



where 7 is the ratio of the constant pressure to the con- 
stant volume specific heats of the gas. In this case wave 
speed scales as the square root of the absolute tempera- 
ture T of the gas and is independent of pressure. For a 
real gas, this is true only in the limit P — > 0. At non-zero 
pressures, a more general theory predicts a dependence of 
c on pressure |32j . Anyway, in the context of the present 
work, no deeper investigation on this subject is needed: 
changes of the value of c will show as a rescaling of the 
frequency spectrum. 



C. Frequency spectrum of a spherical resonator 

In a bounded volume, the homogeneous Helmholtz 
equation associated with Eq. Q determines a discrete 
frequency spectrum. In a spherical geometry the equa- 
tion can be solved by variable separation. The angular 
part of the wave function is given by the spherical har- 
monics Yi m {0,<t>) [33], while the radial function R{r) is a 
solution of 



d 
dr 



,dR 
dr 



= [l(l + l)-k 2 r z ]R 



(6) 



where k = 2m// c — lj/c is the wave vector. Solutions 
of this equation are the spherical Bessel functions of first 
and of second kind 1 33 : 



Ri(r)=Aiji(kr)+B l yi (kr) 



(7) 



with Ai and P>i arbitrary coefficients. 

Let us suppose that the walls of the cavity at radiuses 
r\ and r 2 are perfectly rigid. According to Eq. (J2| , on the 
cavity walls, where the radial component of the velocity 
is zero, the acoustic wave satisfies Neumann boundary 
conditions R[(ri) = J2j(r 2 ) = 0, or 

A l f l (kr 1 )+B l y' l (kr 1 )=0 
A l f l (kr 2 ) + B l y' l (kr 2 )=0. 

If both equations must hold, then j'^kri) y[{kr 2 ) — 
Ui(kri) j'i(kr 2 ) = 0. This last equation generates a dis- 
crete spectrum of allowed wave vectors k n \ , with index n 
numbering the roots for a given I. The solution k = for 
I = is left aside for the purpose of the present work. The 
radial parts of the normal modes of a spherical acoustic 
cavity are then of the form 



Rni(r) = A' ni 



■II, \ fl( k nl r l) s 



We call radial modes the spherically symmetrical ones 
with / = 0, and angular modes the others. 

If 7*1 = 0, namely the cavity is an empty sphere, the 
coefficient P>i in Eq. (]7|) is zero; Neumann boundary con- 
ditions can still be imposed and the set of allowed wave 
vectors changes accordingly 34J. In any case, at a single 
frequency lo 71 i we have 

P(r,u nl ) = P + R nl {r)e-^ t+ ^ A'^ m Y lm (8 , <f>) . 
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Table II. Allowed wave vector values (in inverse meters) of 
acoustic waves inside a gas filled spherical resonator with r2 = 
20 cm, for n = 4 cm and ri = 6.35 cm. 
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Table III. Allowed wave vector values (in inverse meters) of 
acoustic waves inside the gas filled spherical resonator used by 
Dylla and King [22j . having ri — 1.27 cm and T2 = 19.7 cm. 



Tab. [TT] lists the first allowed wave vectors for the two 
cavities used in the present work, both with external ra- 
dius r 2 = 20 cm, and with the internal wall at n = 4 cm 
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and r\ = 6.35 cm, respectively. In Tab. |III| the allowed 
wave vectors of the apparatus used by Dylla and King 
[22] are shown. Note that for this last configuration the 
radial modes are almost degenerate with angular modes, 
what makes that cavity not quite suitable for the experi- 
ment. In fact, the authors meant to tune the (1, 0) mode, 
which is overlapping the (1, 3) for a Q factor w 1000. 



D. Acoustic effect of electric polarization 
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Table IV. Physical parameters used in the calculations. Sound 
speed refers to T = 298 K and Po = 1 atm. Static polariz- 
ability values are taken from Ref. [35] . 
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Table V. Resonance frequencies (in hertz) of the first radial 
modes of two cavities with r 2 = 20 cm and ri = 4 cm and 
ri = 6.35 cm, filled with different gases. All the values refer 
to P = 1 atm and T = 298 K. 



We want now to determine the acoustic effect of elec- 
tric polarization on a gas contained inside a spherical 
capacitor with perfectly rigid walls. In doing this we will 
consider only acoustic modes having spherical symmetry, 
namely radial modes with 1 = 0. The parameters used in 
the calculations are summarized in Tab. IV Tab. |V| lists 
the resonance frequencies of the first radial modes of the 
two cavities used in the present work filled with different 
gases. 



For ideal gases, the divergence of the second term of 
the force density of Eq. ([lj is 
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where 



C = 4ire 



ri - n 



is the capacitance of the spherical capacitor, and V is 
the applied voltage difference. The parameter (3 depends 
on cavity geometry. For cavities with r 2 = 20 cm and 
r*i = 4 cm or fj = 6.35 cm it has values: 
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For harmonic excitations, the wave equation Q re- 
duces to an inhomogeneous scalar Helmholtz equation 
whose solution is obtained by numerical integration with 
boundary condition dP/dr = f r both at r = r\ and at 
r = r 2 ■ In the forcing term given by Eq. ^ the relevant 
quantity is the Fourier component of V 2 at a resonance 
frequency oj n o, which we call V 2 (u> n o). In order to take 
into account the viscosity, we assume a complex wave 
vector: 



(10) 



The relaxation time r depends on base pressure, reso- 
nance mode and cavity; it is related to the quality fac- 
tor through Q = ixv n QT. No attempt has been made to 
deduce r from known values of viscosity, it is instead 
measured independently. 

The output of the calculation is the amplitude of the 
pressure wave P(r, v) at given r and v. The calcula- 
tions are performed for different values of the frequency 
around the resonance frequencies t^„o in order to draw 
the amplitude resonance curve. The results are studied 
as a function of the values of r and of V 2 (oj n o)- From 
the analysis of the results it turns out that the dynamic 
pressure peak values at r 2 due to electric polarization can 
be written as 



P(r 2 ,u n o) = B n0 irv n0 Taf3V 2 (oj n0 ) 
= B n0 Q n0 af3V 2 (uj n0 ). 



(11) 



The results of the calculations for the cavities employed 
in the present work are summarized in Tab. |VI| For 
V ~ 1 kV and Q ~ 10 3 , expected values of the dynamic 
pressure are P ~ 1 mPa. 

A similar calculation performed for the apparatus of 
Ref. [H] (r-r = 1.27 cm, r 2 = 19.7 cm, SF 6 at P w 
0.8 atm, V = 3500 V rms, Q = 1100) gives P{r 2 ,v w ) = 
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Table VI. Numerically calculated results of acoustic effects 
due to electric polarization in spherical cavities having r 2 = 
20 cm, for r\ — 4 cm and n = 6.35 cm. Dynamic pressure 
values are calculated through Eq. |TTJ . All the numbers refer 
to P = 1 atm and T — 298 K - see Eq. jgj. 



1.42 mPa, a value which is about three times larger than 
the 4.68 mdyn/cm 2 that can be obtained from Eq. (3) 
of the same reference. As a consequence, since they use 
the polarization measurements to perform an absolute 
calibration of their instrument, the ratio of the voltage 
output of the microphone to the expected pressure signal, 
V s /Pd, given in Eq. (4) of Ref. 22 , is too large of the 
same factor three. 



E. Acoustic effect of non-neutrality of matter 

The acoustic effect of an hypothetical charge asym- 
metry of proton and electron can be tackled in a similar 
way as the polarization effect. Unlike that case, however, 
for the non neutrality of matter the solution of the wave 



equation (12 1 can be obtained in a closed form. In the 
following we show a simplified version of the analytical 
solution. 

The divergence of the first term of Eq. ([lj is zero ev- 
erywhere inside the capacitor. Hence the wave function 
satisfies the homogeneous wave equation ([6| which, for 
radial modes, is 



1 d 
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r 2 dr \ dr 
with inhomogeneous boundary conditions: 
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where M — p m oRT/ Pq is the molecular mass and we 
have assumed e p _ c = e n = t q . Here, again, the coefficient 
8 does not depend on the gas species, but on the geometry 
of the cavity. Its values for cavities with r-i = 20 cm and 
r*i = 4 cm or r\ = 6.35 cm are: 

e P nrj 8 C>mol _Ffe 298 K 

d 4 = = 1.96 10 ? 

m p RT T2 — T\ kg-m 1 atm T 



(13) 



^6.35 



3.64 10 ; 



a C-mol P 298 K 

o 



kg-m 2 1 atm T 



From Eq. 
P(r)=Aj (kr)+By Q (kr) 



A sin(fcr) B cos(fcr) 



kr 



kr 



The coefficients A and B are determined by the boundary 
conditions 



Aj' Q {kr l ) + By' Q {kr l )=e ll M8 
Aj' (kr 2 ) + By' (kr 2 ) = e q MS 



V 



V 



which give 



a » r n; j cos(fcri) — cos(fcr 2 ) + kr\ smikrA — kr 2 sm(fcr 2 ) 

A = e„ M o V k — — — - — — — — — 

kAr cos(fcAr) — (1 + k l r\r2) sm(fcAr) 

t-> 1 1 r , ^ r i cos(fcri) - kr 2 cos(fcr 2 ) - sin(fcri) + sin(fcr 2 ) 

B = e„ M o V k — — —5 r — , 

— kAr cos(kAr) + (1 + fcVir 2 ) sm(feArj 

where Ar = r 2 — n . In turn one obtains 

P(r) — e M SV cos ^'- r ~ — ^ r2 cos i^( r 2 — r )] + sin[fc(r — n)] + sin[fc(r 2 — r)] 
9 r [fcAr cos(fcAr) — (1 + k 2 rir 2 ) sin(feAr)] 



I 

When the wavevector approaches a resonance value k n o, this expression diverges, as the denominator goes to zero. 
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To take into account dissipation, we use the complex def- 
inition of k derived from Eq. ( 10 1 in the case of weak 
dissipation: 



i 

2Q 



We note that for the real and imaginary parts of the 
wavevector it holds kj <C k^. The complex value of k 
can be substituted into the expression of P(r), whose real 



and imaginary parts are in phase and out of phase with 
respect to the forcing field. To obtain an estimate of the 
amplitude of P{r) at resonance, we consider separately 
the numerator and the denominator of the expression of 
P{r) as a function of k. Let's call them N(k) and D(k), 
respectively. At resonance we have 

N{k n0 + ifcj) N{k n0 ) + ikiN(k n0 ) 2Q N(k n0 ) 



D(k n0 + iki) D{k n0 ) + ikiD(k n0 ) ik n0 D(k nQ ) 
hence 



P{r 2 ,i>na) 
e q Q V(ui n0 ) 



2M5 



k n0 ri cos(fc n0 Ar) - k n0 r 2 + sin(fc n0 Ar) 

2 l0 r ^h k norir 2 Ar cos(k n0 Ar) - (rf + if) sin(fe n0 Ar)] 



(14) 



This result is incompatible (analitically and numerically) 
with Eq. (2) of Ref. [22]. 



n 


value 


species 


mode (n, 0) 


(cm) 






(1,0) (2,0) (3,0) 




D n0 (m" 1 ) 




2.27 


0.557 


0.460 






SF 6 


65.0 


15.9 


13.2 




P(r 2 ) 


Xe 


58.3 


14.3 


11.8 


4.00 


Kr 


37.4 


9.17 


7.57 


e q QV(aj n0 ) 




(MPa/V) 


Ar 


17.8 


4.37 


3.61 






N 2 


12.5 


3.06 


2.52 




D n0 (m" 1 ) 




1.87 


0.311 


0.298 






SF 6 


99.4 


16.5 


15.8 




P(r 2 ) 


Xe 


89.2 


14.8 


14.2 


6.35 


Kr 


57.2 


9.51 


9.11 


e q QV(Lo n0 ) 




(MPa/V) 


Ar 


27.2 


4.53 


4.34 






N 2 


19.1 


3.17 


3.04 



Table VII. Calculated results of acoustic effects due to charge 
asymmetry in spherical cavities having r 2 = 20 cm, for 
ri = 4 cm and r\ = 6.35 cm. Dynamic pressure values are cal- 
culated through Eq. ( 15 1 . All the numbers refer to Po = 1 atm 
and T = 298 K - see Eq.JT3|. 



In analogy with the case of electric polarization, the 
expression of the dynamic pressure peak values at r 2 due 
to the charge asymmetry is cast in the form 

P(r2,Vno) = D n0 Q n0 e q M SV{uj n0 ). (15) 

The results for the caviti es con sidered in the present work 



VII 



If we suppose e q 



lfr 21 , 

1 ^Pa 



are summarized in Tab. 

V ~ 1 kV and Q ~ 10 3 , a dynamic pressure P 
is expected. 

A similar calculation performed on the first radial 
mode of the cavity of Ref. [Hj (ri = 1.27 cm and 



= 19.7 cm, 1 atm of SF 6 
P(r 2 , u w ) 



gives the following result: 

. = 21.2 MPa/V (exp. of Ref. [22] ), 

to be compared with the corresponding value that can 
be calculated using Eq. (2) of the same reference, which 
is a factor about sixty larger (Pq = 1.24 GPa/V). Also 
in this case we suspect that the neglecting by Dylla and 
King of the Bessel functions of the second kind in the 
solution of the wave equation might be the reason for the 
discrepancy. 

III. EXPERIMENTAL SET-UP AND METHOD 

A principle scheme of the experimental apparatus is 
shown in Fig. [2j The core of the set-up is a 40-cm di- 
ameter spherical cavity made of two mechanically joined 
hemispheres carved from a block of A15056 aluminum al- 
loy. The internal surface of the hemispheres is machined 
within a 0.1 mm precision and polished with 6 — 8 /im 
sandpaper. The thickness of the shell is not uniform, as 
roughly sketched in Fig. [3j Two holes are drilled in the 
sphere, a 1.6 cm diameter one located at the north pole 
and a 2 mm one at the south pole. In the course of the 
experiment the two holes played different roles: in a first 
phase, the upper hole was a gas inlet and was crossed by 
the suspension wire of the inner electrode; a bellows oper- 
ated translation stage allowed the vertical positioning of 
the central electrode with a 1 mm resolution; the acous- 
tic signal was recorded at the lower hole. Later, when 
there was evidence that the empty volume of the bellows 
was altering the amplitude of the resonant modes of the 
cavity due to the phenomenon of Helmholtz resonance 
|37) . the height adjusting facility was abandoned and the 
upper hole was closed with an home made Macor electric 
fccdthrough; the lower hole assumed then also the func- 
tion of gas inlet. Two spherical internal electrodes have 
been used in turn: an 8 cm diameter aluminum one sus- 
pended to a 0.25 mm brass thread, and a 12.7 cm stainless 




Figure 2. Principle scheme of the experiment. 



steel one hanging by a 1 mm stainless steel wire. In or- 
der to provide acoustic and thermal isolation, the sphere 
is housed inside a large vacuum chamber where a rough 
vacuum (~ 10~ 2 mbar) is made. Inside the chamber 
the sphere lays on a thick rubber pad. The apparatus 
is supported by a pneumatically damped optical table; 
by varying the relative pressure of the legs, the tilt angle 
could be controlled at the milliradiant level. The cav- 
ity can be evacuated down to 10 -5 mbar and then filled 
with very high purity gases up to a pressure of 10 atm. 
The complete set-up lies in a temperature regulated room 
(T sa 295 K). After gas insertion, a time delay of at least 
one day was generally kept before taking measurements. 

The outer shell is in electric contact with ground and 
isolated from the inner electrode. An HV feedthrough 
connects the inner electrode to an HV amplifier driven 
by a waveform generator; in this way an oscillating elec- 
tric field is established inside the cavity. The HV am- 
plifier is composed of a very high linearity power stage 
capable of up to 100 V pp output (5 V pp input) with an 
output impedance less than 1 f2. The power amplifier is 
followed by a 1:100 transformer wound around an ultra 
linear M0 grade magnetic core, with total power 1.5 kVA. 
The whole HV amplifier is actually located in a different 
room than the sphere, inside a Faraday chamber. The 
ratio of the amplitudes of the second to the first har- 
monic has been measured to be less than 5 x 10~ 5 at the 
amplifier output. Ac and dc components of the electric 
voltage Vhv are measured using two six and a half digit 
multimeters. The ac readout is made through a 1:1000 
voltage divider. 




Figure 3. Top: schematic mechanical drawing of the cavity; 
the two different values of the internal electrode diameter are 
drawn. Bottom: detail of the microphone mounting. 



As the microphone is sensitive to electric fields, it could 
not be flush mounted, namely with the sensitive mem- 
brane mimicking the inner surface of the cavity. In fact, 
we observed that in this last configuration an oscillating 
electric field generates a spurious acoustic signal which 
grows with the square of E. Acoustic oscillations of the 
gas inside the cavity are measured instead through a 
2 mm diameter, 5 mm deep hole by means of an half 
inch prepolarized preamplified free-field condenser mi- 
crophone (PCB 378B02); the head of the microphone 
is plugged inside a Teflon socket which lays against the 
outer surface of the sphere; an half-inch diameter Viton 
o-ring seals the contact between Teflon and aluminum 
(see lower part of Fig. [3]) ; the o-rings are not meant to 
stand any differential pressure, as the whole microphone 
assembly is in the same atmosphere of the cavity. The 
microphone is mounted on a linear translation stage; this 
measure became necessary since when the south pole hole 
was the only opening left in the sphere, to allow for vac- 
uum pumping and gas filling the cavity. 

Nominal microphone sensitivity is K m = 50 mV/Pa at 
1 atm in the frequency interval 3.15 Hz - 20 kHz; each 
microphone employed in the experiment had a calibra- 
tion sheet, the actual value of the sensitivity being known 
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Figure 4. (Color online). Measured ratio Gh of the fre- 
quency response in air of the microphone mounted as in 
Fig. [3j to its free field response for three different distances 
from the loudspeaker: 25 cm, 35 cm and 58 cm. As the 
three graphs are quite similar, the smooth curve fits the 
average (not shown) of the three with the resonance func- 
tion y = A[{v 2 - vl) 2 + v 2 T 2 }~ l/2 with A = 4.58 kHz 2 , 
uq — 2.13 kHz, r = 0.80 kHz. The arrows mark the positions 
of the first radial modes (n, 0) for the cavities with n = 4 cm 
and ri = 6.35 cm. 



with a 0.01 mV/Pa uncertainty. At higher (lower) pres- 
sures the microphone loses (gains) sensitivity at a rate of 
0.013 dB/kPa, so that K rn = 27.3 mV/Pa at 5 atm. The 
noise limit of the preamplifier is 1 /iPa/Vflz. The cavity 
and the duct in front of the microphone constitute an 
Helmoltz resonator, whose response has been measured 
in air and compared to the free held response of the mi- 
crophone. Fig. |4]shows the results of this calibration: the 
Helmholtz cavity amplifies a factor Gh ~ 1-6 at the first 
radial mode of the cavity with r\ = 6.35 cm, much the 
same as the second radial mode, while the third radial 
mode is attenuated. The microphone is also sensitive to 
vibrations, hence it can be used to record a spectrum of 
the normal modes of vibration of the aluminum skin. In 
Fig. [5] such a spectrum is shown; the data are obtained 
by mechanically exciting the evacuated aluminum cavity 
with an external source. As can be seen, no mechanical 
resonance appears below 1 kHz. 

The microphone output is amplified and read using 
a dual phase Lock-In Amplifier (LIA) referenced to the 
waveform generator. A FFT spectrum analyzer is also 
used to analyze the microphone signal. To avoid electric 
pick-ups the LIA reference input is connected to the sync 
output of the signal generator through an optical link. 
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Figure 5. Acoustic signal from the microphone due to me- 
chanical vibrations of the evacuated cavity. A small loud- 
speaker is used as an excitation source. 



This system allows pressure measurements of the type: 

P(t) = Pi cos(uit + <j>), (16) 

where the dynamic pressure amplitude Pi and its phase 4> 
are directly calculated from the (amplitude, phase) out- 
puts of the LIA or of the FFT analyzer. 
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Figure 6. Decay of the demodulated microphone signal from 
the (1,0) resonance mode of the cavity with inner electrode 
radius n = 6.35 cm filled with 1 atm of SF6; the exponen- 
tially decaying fit line, completely masked by the experimen- 
tal points, gives a quality factor Q — ■KVxar = 4780. 



In order to exploit the amplification factor of the acous- 
tic modes, all the measurements have been performed 
in correspondence of the first cavity radial modes, hav- 
ing resonance frequencies v n Q. The acoustic modes of 
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Figure 7. (Color online). Lineshape and phase of the (1,0) 
resonance mode of the cavity with inner electrode radius 
n = 6.35 cm filled with 1 atm of SFg; the two lines are in- 
dependently fit, giving compatible Q values; the fit functions 
are completely masked by the experimental points. 




2 3 4 

Pressure (bar) 

Figure 8. Pressure dependence of the quality factor of the 
(1,0) resonance mode for SFg in the cavity with ri = 6.35 cm 
(upper curve, diamonds) and for Ar in the cavity with rj = 
4 cm cavity (lower curve, squares). Continuous lines are only 
meant to guide the eye. 



the cavity are probed exciting sound waves either using 
a loudspeaker external to the sphere or through electric 
polarization, by placing an oscillating high voltage on the 
inner electrode. The frequency of the resonant modes is 
determined maximizing the output of the LIA and mon- 
itored through its phase. The associated quality factor 
Q is obtained from the decay time of the LIA output 
when the exciting source is switched off (see Fig. [6]) . By 
scanning the exciting frequency, the lineshape of the res- 
onance mode can also be traced directly from the output 
of the LIA or of the FFT (Fig. [7}. The quality factor 
of a mode depends on pressure, increasing with increas- 
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Figure 9. Observed spectrum and Q values of the cavity with 
n = 4 cm, with 1.2 bar of SFs. 
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Figure 10. Resonance frequency of the (1,0) mode of the 
n = 4 cm cavity filled with 1 atm of N2 as a function of the 
vertical position of the central electrode. Following Ref. [31] 
the data are fitted with a parabolic curve. 



ing pressure (Fig. ^| . The measured frequency spectrum 
of the cavity (Fig.T9l) shows an excellent agreement with 
the calculated spectrum. Frequency shifts due to finite 
admittance of the cavity walls or to geometrical imper- 
fections |38j are beyond the experimental sensitivity and 
the scope of the present work. The frequencies of the 
acoustic modes are seen to depend on the position of 
the central electrode. In Fig. 10 we have an example 
of this behavior: the resonance frequency shows a maxi- 



mum which corresponds to a minimum in the eccentricity 
of the electrode 1531. 



Measurement type 


VHV 


V 


Vi 


P phase <j) 


V 2 (w n0 ) 


Polarizability I (PI) 


VnO 


^0 


^0 


4>HV — 7r/2 


2V0V1 


Polarizability II (PII) 


V n0 /2 







2(f>HV — 7r/2 


Vf/2 


Neutrality (N) 


V„o 





¥=0 


(j>HV — 7r/2 





Table VIII. Experimental parameters for the different con- 
figurations of measure. The Fourier component of V 2 at a 
resonance frequency uj„o, V 2 (uj n o), does not apply to neutral- 
ity measurements, which depend linearly on V. 
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The measurements of interest are those connected with 
the presence of the electric field in the resonator. The 
high voltage signal on the inner electrode is parametrized 
as: 

Vffv(t) = V + Vi e - i ( u *vt+*Hv) < 

Three types of measurement are possible, depending on 
the values of Vo and V\ ; all the relevant values of the ex- 
perimental parameters are summarized in Tab. |VliT} Like 
in the original experiment by Dylla and King, the po- 
larizability measurements are used to test and calibrate 
the apparatus. When the high voltage excitation fre- 
quency coincides with the frequency of a resonant mode 
and Vo ^ 0, both effects due to polarization and to an 
hypothetical non neutrality of matter are present. In 
this case, considering the pressure signal as a function of 
Vo, one sees that the neutrality measurement N can be 
obtained as the limit for Vo — > of the polarization mea- 
surement PI. Moreover, this measurement configuration 
allows also a continuous monitoring of the forcing field, 
thus reducing phase drift effects. 

IV. RESULTS AND DISCUSSION 




Figure 12. PH-type polarization measurements for several 
SF6 pressures on the (1,0) resonance mode of the cavity with 
n = 6.35 cm, as a function of V\. The data are fitted with Vi 
functions. Resonance frequencies and quality factors range 
from via = 550 Hz and Q = 1400 at P = 260 mbar to 
vio = 537 Hz and Q = 9800 at Po = 3.01 bar. Uncertainties 
between 20 and 50 nV result in reduced x 2 ' s equal to one. 



A. Polarizability measurements 



- 




V (V) 



Figure 11. Typical Pi-type polarization measurement for 
0.25 bar of SFg on the (1,0) resonance mode at v — 551.6 Hz 
of the cavity with n = 6.35 cm, for Vi = 785 V rms, as a 
function of Vo. All the 13 points drawn lay on a continuous 
straight line having slope 32.3 x 10~ 9 and intercept 49.8 nV. 
The quality factor of the cavity is Q = 1400. A 25 nV un- 
certainty on all the points results in a reduced x 2 equal to 



In Figs. [TT| and [T2| typical results of the polarizability 
measurements are shown. In both cases no deviation of 
the data from the theoretical fit functions is observed. 
The slopes and curvatures obtained from the fits allow 
to calculate the electric polarizability of the gases. The 
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Figure 13. Measurements of the molecular polarization a 
of SFg as a function of pressure compared to the value 
a C RC = 7.277 10" 40 C-m 2 /V. In this graph the 0.013 dB/kPa 
correction of the sensitivity of the microphone has not been 
taken into account. 



data scale well with the gas pressure and with the reso- 
nance mode. A normalized plot of determinations of a in 
SF6 at different pressures is shown in Fig. [T3j the mea- 
surements have been obtained from linear fits (Pl-type), 
quadratic fits (PH-type), and from single measurements 
taken in a period of several months during which the ap- 
paratus underwent major changes that affected the un- 
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certainty of the measurements but apparently did not 
alter the data distribution. As can be seen, the exper- 
imental value of a is a factor R a sw 1.25 larger than 
the tabulated one. The reasons for that arc not clear. 
Note however that, since the polarization measurements 
represent a calibration of the apparatus at the working 
frequency, we have assumed that all the acoustic signals 
detected by the apparatus, including those measured in 
neutrality measurements, are sensed with a microphone 
sensitivity K' m = R a GHK m 100 mV/Pa. In the fol- 
lowing we consider two possible explanations for the dis- 
crepancy. 

First of all, the gas is not isolated from the aluminum 
skin of the cavity. In fact, as mentioned before, a loud- 
speaker external to the cavity is able to excite the reso- 
nance modes of the gas. Conversely, during P-type mea- 
surements, an accelerometer in external contact with the 
cavity is able to detect the vibration of the wall of the 
cavity and decay curves of the type of Fig. [6] can even 
be recorded; no such effect can be observed if vhv is 
detuned from resonance, meaning that the coupling be- 
tween the high voltage and the cavity skin is mediated 
by the resonating gas: the gas-skin system should then 
be considered a single resonating system. In the origi- 
nal experiment by Dylla and King, the cavity skin was a 
thin foil of copper, a less rigid material. A second possible 
source of the difference between calculation and measure- 
ment is the fact that the presence of the electric wire, to 
which the inner electrode is suspended, makes the elec- 
tric field lines not perfectly radial; hence the acoustic 
and the electric fields do not completely match. More- 
over, the highest values of the electric field are found near 
the wire. 

We also performed measurements of the polarizability 
of other gases at 1 atm. For Xe and Kr the relative 
polarizability a re i = a/«cRC is similar to the value for 
SFg: a re i = 1.29 for Kr and a rc i = 1.25 for Xe. As the 
resonance frequency increases, we find a rc i = 1.7 for Ar 
and a rc i = 1.55 for N2, due to the onset of the mechanical 
resonances (see Fig. [5]). 



B. Neutrality measurements 

As far as the neutrality measurements are concerned, 
no direct N-type measurement was possible. In fact, 
when Vq = and vhv = ^n0> a small but significant 
spurious signal is always present, with the right phase 
and an amplitude corresponding to a non-neutrality of 
the matter of the order of e q ~ 10 -19 . This signal is not 
an electric pick-up, as is proved for example by detuning 
the excitation frequency; moreover, it decays with the 
time constant of the cavity. It can have a mechanical 
contribution, due to waves excited by the high voltage in 
the skin of the cavity, an effect completely negligible as 
long as polarizability effects are concerned, but that can 
play a role in this context. A large contribution must 
however come from free charges in the gas. In fact, we 
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Figure 14. Demodulated dynamic pressure signal (y or in 
quadrature output of the LIA), recorded in N-type measure- 
ment configuration with 0.48 bar of SF6, following the appli- 
cation to the inner electrode of a dc high voltage Vo — 350 V 
for ten minutes. After switching off the high voltage at t = 0, 
Vo decays with a time constant of about one minute. The x 
or in phase output of the LIA is not shown, as during the 
observation time this last signal is fairly stable with values 
near zero. 



observed an increase of the neutrality signal by irradi- 
ateing the cavity with 7-rays from a radioactive source. 
The spurious signal reacts with hysteresis to dc high volt- 
ages applied to the inner electrode for time periods of the 
order of a few minutes. A typical example of this behav- 
ior is shown in Fig. 
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We have also observed that the 
spurious signal grows more than linearly with the base 
pressure inside the resonator. The pressure inside the 
resonator should then be kept as low as possible. With 
this finding, the idea of increasing the pressure inside the 
cavity to obtain larger signals reveals useless. None of 
these effects are reported in the 1973 paper by Dylla and 
King. 

To go around the impossibility of a direct N-type 
measurement, we chose to use instead the intercepts in 
Pl-type measurements taken at a base pressure Pq sa 
0.25 bar of SF 6 . As noted before, the intercepts repre- 
sent a determination of the acoustic effect of an hypo- 
thetical charge asymmetry. We expected that the con- 
tinuous presence of a DC voltage could sweep away the 
free charges from the gas. The hypothesis proved to be 
right: referring to Fig. |11| where a typical set of mea- 
surements is shown, one can see that the value of the 
intercept is much smaller than that obtained in a direct 
N-type measurement. All the measurements have been 
taken on the (1,0) resonant mode at Vw ~ 552 Hz, with 
Q = 1400. The data taking procedure is as follow: the 
slope and intercept of the straight line have been deter- 
mined several times in a period of several weeks; each 
data set has about ten points with at least five different 
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voltage values; the points have been sampled in a non 
monotone sequence with a few minutes between the two 
of them, in order to avoid transient phenomena; the posi- 
tive and negative branches of each line have been sampled 
in different data sets, to keep constant the polarity of the 
dc voltage during each measurement. For each point the 
resonance frequency is tuned observing the phase of the 
polarization signal; the maximum error on the tuning 
of the resonance frequency is less than 4 mHz; for the 
0.25 bar cavity this implies an uncertainty on the ampli- 
tude of 0.04%. A linear regression was performed on each 
data set; all the points of each set were given the same 
statistical uncertainty chosen in such a way to have a re- 
duced x 2 of the order of unity. The calculated statistical 
uncertainties lay within a factor two from the measured 
noise figure of the electronic chain. 
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Figure 15. (Color online). Plot of thirtysix determinations 
of the charge asymmetry e q obtained from the intercepts of 
Pi-type measurements taken on the resonance mode (1,0) of 
the cavity with rj = 6.35 cm with Po — 0.25 bar of SF6. 
The factor R a Gn ~ 2 (see Figs. [4] and 13 1 has been taken 
into account. The continuos lines are the averages, taken 
separately, of the data measured with negative (upper hor- 
izontal line, squares) and positive (lower horizontal line, di- 
amonds) Vo. The two average values are (+9.1 ± 1.4) 10~ 21 
and (—9.3 ± 1.6) 10 -21 , respectively. 

All the results obtained for the value of the intercept 
are shown in Fig. |15| In Fig. [16] we show the histograms 
of the results. We note that not all the values are com- 
patible with zero. It is observed that all but a few data 
sets taken with Vq < have positive intercept, while the 
opposite is true for the data sets taken with Vo > 0. 
If the two groups of intercepts are separately averaged, 
opposite values are found: this means that besides the 
effect of polarization, at least another spurious effect of 
unknown origin must be present which, at variance with 
the hypothetical non neutrality of the matter, changes 



Charge asymmetry xlO 



Figure 16. (Color online). Histograms of the two groups of 
data of Fig. [L5] 



sign with the sign of the excitation signal. We get rid 
of this unwanted effect, thus extracting the value of the 
charge asymmetry, by making the half-sum of the two val- 
ues. Finally, we find for the charge asymmetry of electron 
and proton as well as for the charge of the neutron 



-0.1 ±1.1) 10 



-21 



This result is compatible with zero and the sensitivity is 
at the level of the best results of Tab. Q] 

Improving this result depends little on having better 
statistics. Several other effects concur instead in frustrat- 
ing further work on this set-up: the temperature stability 
of the whole apparatus, the temperature sensitivity of the 
voltage divider used in the readout of the ac component 
of Vhv, the sensitivity and linearity of the multimeter 
employed, the sensitivity of the microphone, the dynamic 
range of the spectrum analyzer. A better understanding 
of the spurious effects observed would also be necessary. 



V. CONCLUSIONS 

The acoustic technique of the 1973 experiment by 
Dylla and King [55] has been employed to test the neu- 
trality of matter. A re-analysis of Ref. [55J shows that 
their claimed limit is invalidated by crucial mathemati- 
cal errors. In our measurement no evidence of a charge 
asymmetry has been found. Assuming charge conserva- 
tion in the /3 decay, our result is at the level of the best 
existing limits, namely 



e,<lx 10~ 21 . 



In the framework adopted, this limit holds both for 
the charge difference of proton and electron and for the 
charge of the neutron. 
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